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Abstract 

The structure of small mixed helium clusters doped with one calcium atom has been determined 
within the diffusion Monte Carlo framework. The results show that the calcium atom sits at the 
4 He- 3 He interface. This is in agreement with previous studies, both experimental and theoretical, 
performed for large clusters. A comparison between the results obtained for the largest cluster we 
have considered for each isotope shows a clear tendency of the Ca atom to reside in a deep dimple 
at the surface of the cluster for 4 He clusters, and to become fully solvated for 3 He clusters. We 
have calculated the absorption spectrum of Ca around the AsAp <— 4s 2 transition and have found 
that it is blue-shifted from that of the free-atom transition by an amount that depends on the size 
and composition of the cluster. 

PACS numbers: 36.40.-c, 33.20.Kf, 67.60.gj 
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I. INTRODUCTION 



Helium clusters are weakly bound quantum systems as a consequence of the small atomic 
mass and the weak van der Waals interaction between helium atoms. There are several recent 
papers [l, 2, 3, Q, Q, H] reviewing the most relevant aspects of the physics and chemistry of 
these remarkable systems, to which we refer the interested reader. 

It was found about twenty years ago that helium clusters could pick up closed shell 
molecules singly and that these molecules were located in the interior of the clusters 0, Q] ■ 
It is now well established that helium clusters may virtually capture any kind of atoms 
and molecules, most of which are located in their bulk. Alkali atoms are the sole known 
exception, as they remain on the surface of the cluster upon their capture. Electronic 
spectroscopy allows to obtain information on the cluster structure, and in particular on the 
position of embedded impurities in the cluster, via the shift and width of the electronic 
transitions 0, 9]. 

Heavy alkaline-earth atoms but magnesium exhibit a different behavior for each helium 
isotope, in contrast with most impurities, which solvate similarly in both. They are solvated 
in 3 He but not in 4 He clusters, in which case they reside in dimples on the cluster surface. 
Although still debated, the magnesium atom is most likely in a very delocalized state within 
the bulk of the 4 He cluster. 

Clusters formed by a mixture of both helium isotopes are particularly interesting, as they 
are made of bosons and fermions with different mass interacting through the same potential. 
It is known that in the bulk, the excess in kinetic energy pushes the 3 He atoms to the liquid 
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It has been 



surface [10|], an effect that also shows up in mixed droplets, see e.g. Ref. 
experimentally found that, depending on the size and composition of the drop, Ca atoms 
reside in bubbles at or near the 4 He- 3 He interface This scenario has been confirmed by 
density functional (DF) calculations Ij), [l^] . Heavy alkaline-earth atoms thus may offer a 
unique opportunity to study the interface of isotopically mixed helium clusters |l4j. 

In this paper we address the study of small mixed 4 He- 3 He clusters doped with one Ca 
atom. Our calculations are based on diffusion Monte Carlo (DMC) methods using current 
realistic interatomic interactions. In particular, we have calculated density distribution 
functions and have estimated the atomic shifts and widths for Ca attached to a cluster, for 
several combinations of the number of isotopes in it. Results for some 4 He and 3 He clusters 
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doped with Ca are also presented to complete the discussion. It is worth stressing that, 
while 4 He clusters doped with different impurities have been thoroughly studied, theoretical 
works on doped 3 He clusters are scarce [la . Ilq . 1 171 I la . Il9| apart from those carried out 



within the DF approach. The paper is organized as follows. In Sec. HI] we give some details 
about the DMC calculations. In Sec. Hill we present our results, and in Sec. IIVI we present 
a brief summary. 



II. METHOD 

The DMC description is based on a variational or importance sampling wave function. 
We have used a rather simple form, which contains the basic required properties. It is a 
generalization of the trial function adopted in our previous studies on pristine mixed helium 
clusters 
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21]. It is written product of seven terms 



(1) 



consisting of a Jastrow factor for each pair of different atoms, plus the spin-up and -down 
Slater determinants required to satisfy the Pauli exclusion principle for 3 He fermions. Here 
{7Z} represents the set of 3 (A/4 + A^ 3 + 1) coordinates of the atoms forming the cluster. Any 
of the Jastrow terms has the generic form 

— atMN r ij 



MN 





i>MN 







(2) 



where indices i,j run over the corresponding type of atom, and includes a short-range repul- 
sion term with parameters bMN and vmn associated to it, and a long-range confining term 
with corresponding parameter olmn- The so defined ^ mn function is explicitly symmetric 
under the exchange of particles. 

The antisymmetry required for 3 He fermions is incorporated in the Slater determinants 
and related to the spin- up and -down fermions. These Slater determinants are of 
primary relevance because they define the set of spin-up and -down nodal surfaces which 
strongly constrain the DMC algorithm. As in our previous works 20j, |21] we have assumed a 
shell-model like structure, taking the single-particle orbitals as harmonic polynomials of the 
fermionic cartesian coordinates, thus guaranteeing that the resulting Slater determinants 
are translationally invariant. Moreover, we have always assumed a filling scheme in which 



the total spin is minimum, either or 1/2, respectively, for N3 even or odd. Also, the so 



called Feynman-Cohen back-flow 



221 ] has been incorporated into the scheme by substituting 



i i = r i + ^2r]( y r ij )(r i -r j 

#3 



(3) 



in the Slater determinants 
form used in Ref. 
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231 ] . For the backflow function 77 (r) we choose the medium-range 



namely 77 (r) = A/r 3 , with the same value of A = 5 A 3 . 



The form of the short-range repulsion term was introduced long ago by McMillan 25| for 
describing the homogeneous liquid 4 He using a 12-6 Lennard- Jones interaction. The values 
of the two parameters v and b are fixed by the short-range behavior of a pair of atoms. In 



our calcu 
Hinde 



potential, and that of 



13. For all cluster sizes 



ations we have employed the He- He Aziz HFD-B(HE) 
271 ] for the X X E Ca-He interaction, as parametrized in Ref. 
we have used the following values: vmn = 5.2 A, 6 44 = 2.95 A, 6 33 = 2.85 A, 5 34 = 2.90 A, 
as in previous works on pristine and isotopically mixed helium clusters, and 6 4 c a — 3.25 A, 
and &3Ca = 3.00 A. Thus, the trial or importance sampling wave function contains only 
five free parameters, namely a 44 , 033, «3 4 , a 4 c a , and a^ca, which have been determined by 
minimizing the expectation value of the Hamiltonian. 



For the DMC algorithm we have used the short-time Green function approximation 28, 
291 ] with an 0(r 3 ) form 30] . It is worth mentioning that this approximate Green function 
satisfies the microreversibility condition (detailed balance condition). Moreover, the random 
process was constrained so as not to traverse the nodal surfaces, using the so called fixed 
node approximation. 



III. RESULTS 



A. Structure and energetics of Ca@ 4 HeAr 4 -|- 3 HeAr 3 clusters 

In Figures [QH] are plotted the densities in the y = plane of the (A 4 ,A 3 ) = (20,20), 
(20,8), (8,20), and (8,8) clusters. After a simulation running for a long thermalization time 
we have stored a large number of walkers (typically 10 6 ) to get a mixed estimator of the 
probability densities. A series of angular rotations have been performed on each walker to 
place the Ca atom on the z-axis, defined as the line joining the impurity and the cluster 
center-of-mass. Afterwards, a projection on the y = plane has been performed. The 
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resulting plots thus exhibit axial symmetry. 

These clusters display an already well known structure, with a core of 4 He atoms sur- 
rounded by a shell of 3 He atoms. Interestingly, in spite of the small number of helium atoms, 
the tendency of Ca to reside at the 4 He- 3 He interface is clearly visible. Note also the peak 
in the 4 He densities near to Ca, and those appearing in the 3 He shell: the one near the 
impurity is reminiscent of the solvation shell that would fully develop in a larger cluster 
made of either kind of atoms. We attribute the 3 He peak distant from the impurity to the 



tendency of the 3 He atoms to reside, whenever possible, far from the impurity (31( because 



the Ca-He interaction is weaker than the He-He one. Within DF theory, a qualitati vely 



similar structure has been found for clusters with A^ 4 = 50 and N 3 = 18, 32, 50, and 68 [13 ]. 

The calculated ground state energies of the clusters, as well as the solvation energy of 
the dopant, defined as 

S(Ca) = £(Ca@ 4 H eA r 4 + 3 He^) - E^Yie Ni + 3 He^) , (4) 



are given in Table [H We mention that Elhiyani and Lewerenz 32j have carried out calcula- 
tions for Ca@ 4 He7v 4 clusters using the same He-He interaction than ours, and a Ca-He pair 
potential they have obtained at the CCSD(T) level of accuracy, as Hinde's one [2^]. Their 
Ca solvation energies for clusters with N4 = 20 and 40 are -12.29 (1) and -16.55 (6) K, 
respectively, in good agreement with ours as displayed in Table HI 

In Fig. [5] are plotted the Ca solvation energies per helium atom as a function of the total 
number of helium atoms in the cluster. Note the conspicuous shell oscillations for N4 = 0. 
We discuss them below. 



B. The case of Ca@ 3 Hejy 3 clusters 



The study of the lightest 3 He clusters has a great interest, as they are a challenge for 
theoretical methods due to their very small binding energy In particular, pure 3 He clusters 
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and small mixed 4 He- 3 He 
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361 ] . It has been shown 



are bound only for sizes greater than N3 ~ 30 atoms 
clusters present instability regions as a function of N± and N3 [2JJ, 
that the presence of one He atom is enough to trigger the binding of twenty 3 He atoms [21], 
showing the importance of the zero-point motion in the stability of small helium clusters. 
The Ca-He interaction potential is weaker than the He-He one, and its minimum lies at 
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a larger distance. Nevertheless, due to the different zero-point energy of Ca and He atoms, 
we have found that the presence of one single Ca atom causes that clusters made of any 
number of atoms are bound, irrespective of their bosonic or fermionic character. This is an 
obvious result for 4 He clusters, but it is not for 3 He clusters. We give in Table [TT] the total 
energy of Ca@ 3 He7v 3 clusters, and have already represented them, divided by N 3 , in Fig. \5\ 
(N4 = results) because as the pristine 3 He clusters are not bound, it makes sense to plot 
these energies together with the solvation energy of Ca discussed before. It is interesting 
to see the sharp minima at N 3 = 2, 8, and 20. These numbers correspond to the first 
shell closures of the three-dimensional spherical harmonic oscillator, in spite that the helium 
density distributions are far from being spherical, especially for small N 3 values (see Fig. [7j). 
We have no plausible exp lanation for this finding. However, it is worth recalling that in 
cylindrical symmetry 371 ]. the first two minima also arise from the la and la, 2cr, In shell 



closures corresponding to this geometry, see e.g., Ref. [13. 

For Ca@ 3 He7v 3 clusters we may also calculate the differences 



A(JV 3 ) = E(Ca@ 3 Re N3 ) - E(Ca@ 3 He 7V3 _ 1 ) (5) 



than can be interpreted as the negative of the He evaporation energies [38|. We have plotted 
—A(N 3 ) in Fig. EJ The large drops in the evaporation energy after N 3 = 2, 8, and 20 correlate 
well with the "magicity" of these values, i.e., with shell closures. In the N 3 = 9 — 20 range, 
the evaporation energy displays some structure, with a weak local maximum at N 3 = 12. 

As representative examples, the density distributions in the y = plane for the Ca@ 3 Hejv 3 
clusters with iV 3 = 2, 3, 4, 5, 8, and 20 are plotted in Fig. [71 It can be clearly seen the tendency 
of the Ca atom to become solvated in 3 He clusters. It is also interesting to compare the 
morphology of 4 He and 3 He clusters doped with Ca. To this end, we have plotted in Fig. [8] 
the density distributions corresponding to helium clusters made of 40 atoms doped with it. 
It can be seen that Ca is in a dimple state on the surface of the 4 He cluster, whereas it 
moves towards the bulk of the cluster in the case of 3 He. In spite of this bulk location, the 
Ca@ 3 He 40 cluster is not spherically symmetric. This also happens for clusters doped with 
other weakly attractive impurities, such as Mg in small 4 He clusters [38| . The reason is the 
tendency of the system to have a helium density as close as possible to its liquid value in 
order to maximize its binding. This is easier to archieve in an off-center location for small 
clusters, even for N 3 = 2. In this case the two 3 He atoms are slightly apart as a consequence 
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of their large zero-point motion, not fully compensated by the attraction of the Ca impurity. 
This is the reason of the very low He density in a cylindrical region surrounding the z-axis, 
see the corresponding panel in Fig. [3, which is absent in the other clusters. Eventually, 
when the number of 3 He atoms increases, the center-of-masses of the helium moiety and the 
impurity nearly coincide 



C. Absorption spectrum of Ca in 4 Hejv 4 + 3 He7v 3 clusters 



The DMC calculation provides us with a set {71} of walkers indicating the instantaneous 
position of each atom in the cluster. From these walkers we have calculated the density 
distribution functions plotted in the previous subsections. Here we use them to obtain the 
absorption spectrum of the 4s4p l P\ <— 4s 2 1 S$ transition in Ca. 

Electronic spectroscopy is a powerful tool to disclose the structure of impurities in helium 
clusters, since the shift and width of the electronic transitions is very sensitive to the dopant 

n n n 

environment [3|, |5j. Lax method [39| has been since long time ago the standard way to 
determine the absorption spectrum of a dopant atom in helium clusters. It makes use of the 
Franck-Condon principle within a semiclassical approach. It has been adapted and used by 
several authors to analyze the absorption spectrum of e.g., lithium in solid H 2 |4(jj, and of 



different atoms in helium clusters 
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42 



43 



45| . The particular case of Ca atoms we are 



interested in has been addressed within DF in pure and mixed helium clusters 12|, |43J. Lax 



46|, in 



method is usually applied in conjunction with the diatomics- in- molecules approach 
which the atom-cluster complex is treated as a diatomic molecule, the helium moiety playing 
the role of the other atom. 

We have calculate the line shape of the electronic absorptium transition in Ca as 



/Hex J dn\^ gs (TZ)\ 2 5(uj + V gs (TZ)-V ex (TZ)) , 



(6) 



where {TZ} refers to the positions of the atoms, and V gs and V ex are, respectively, the ground 



and excited states potential energy surfaces. We have used the Ca- 



rt cf. 



271 for the ground state, and the 1 U and X S potentials of Ref. 
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3e X X E interaction of 
for the excited states. 



In short, for a given value of u, a 3 x 3 matrix has to be diagonalized to determine the 
three components of the absorption line, each one arising from a different potential energy 



surface, i.e., eigenvalue of the excited energy matrix, see Refs. 



44 



for the details. 
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In Fig. [9] is plotted the absorption spectrum of Ca for several combinations of N± and A3 
values. The energies are referred to the free Ca atom value (23650 cm -1 ). For so few helium 
atoms, there is no appreciable shift, but an appreciable width increasing as A3 does at fixed 
A 4 . A large enough number of walkers bear the information corresponding to a quasi- free 
Ca atom within the simulation volume. This is the reason of the narrow peak at u = 0. 
A similar peak appears in the experiments, and it is due to the excitation of gas-phase Ca 



atoms present in the doping chamber 



47|]. 



In Fig. [10] is plotted the absorption spectrum of Ca for three selected (JV4, A 3 ) combi- 
nations with A4 + A3 = 40. The contributions from the three potential energy surfaces 
are shown, and the long tail at high frequencies is identified as likely arising from the very 
repulsive contribution of the X S pair potential. The shift in the (0, 40) case is clearly seen. 
Finally, in Fig. [TTJ we display the absorption spectrum of Ca@ 3 HeAr 3 for increasing values of 
A 3 . Not surprisingly, the A 3 = 2 case shows no appreciable shift and a small width. As A 3 
increases, both the shift and width progressively do. 

IV. SUMMARY 

Within the diffusion Monte Carlo framework, we have studied the structure and energetics 
of small mixed 4 He- 3 He clusters doped with Ca. We have found that a single Ca atom is 
able to produce bound clusters made of any number of helium atoms, irrespective of their 
fermionic or bosonic character. In the case of 3 He clusters, this is a somewhat unexpected 
result, as the Ca-He interaction is weaker than the He- He one. We attribute it to a zero-point 
motion effect. 

We have found that Ca resides in a deep dimple at the surface of 4 He clusters, and 
in the bulk of 3 He clusters, in agreement with experiments carried out for large helium 
clusters 

In the case of mixed clusters, our calculations show the tendency of the impurity to sit at 
the 4 He- 3 He interface, although the small number of fermions that we can microscopically 
handle does not allow us to make a clear distinction between the interface and the rest 
of the mixed cluster^ Morphologically, we have found that doped mixed clusters obtained 
within DF theory [3] are qualitatively similar to those obtained within the microscopic 
DMC approach. 



S 



Finally, we have used the DMC walkers to semiclassically obtain the absorption spectrum 
of Ca around the 4s4p <— 4s 2 transition, and have found that it is slightly blue-shifted from 
that of free-atom transition, with a clear dependence on the size and composition of the 
mixed cluster in spite of the small number of helium atoms in the studied systems. We 
want to mention that the excitation and emission spectra of Ca atoms implanted in liquid 
4 He and 3 He have been measured [48|] , and that their determination in the case of liquid 
mixtures could also be carried out as a function of pressure and composition. We are at 



present extending the method implemented in Ref. 



12 1 to address this issue. 
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TABLE I: Ground state energies (in units of K) of mixed helium clusters, pristine and Ca doped, 
and solvation energy of Ca for several combinations of N4 and 7V 3 . 
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TABLE II: Ground state energies of Ca@ 3 Hejv 3 clusters 
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FIG. 1: Density distributions in the y = plane for the Ca@ 4 He2o+ 3 He2o cluster. Left: 3D 
representation. Right: grey scale plot of 3 He (top) and 4 He (bottom) densities. The darker the 
region, the higher the helium density. 
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FIG. 5: (color online) Solvation energy of Ca per helium atom for several Ca@ 4 HeAr 4 + 3 He7v 3 
clusters, and total energy per helium atom of Ca@ 3 HeAr 3 as a function of the total number of 
helium atoms in the cluster. 
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FIG. 6: Evaporation energies of 3 He atoms in Ca@ 3 He7v 3 clusters for the indicated N3 values. 
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FIG. 8: Grey scale plot of the density distributions in the y = plane for Ca@ 4 He4o and Ca@ 3 He4Q. 
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FIG. 9: (color online) Absorption spectrum of Ca in mixed clusters with N4 = 8 (left panel) and 
-/V4 = 20 (right panel), and three N$ values. 
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FIG. 10: (color online) Absorption spectrum of Ca for three selected combinations of helium atoms 
in clusters with ./V4 + ./V3 = 40. 
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FIG. 11: (color online) Absorption spectrum of Ca in 3 Hejv 3 clusters for the indicated number of 
atoms. For clarity, the N% = 2 and 40 results have been drawn in a separate panel. 
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